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Report  on  Final  Results 

The  Final  Report  is  based  around  the  central  task  put  forth  in  the  original  proposal:  to 
theoretically  estimate  and  optimize  the  figure  of  merit  ZT  for  the  type  II  superlattice  (T2SL), 
and  estimate  whether  it  might  work  as  a  cryogenic  Peltier  cooler.  This  theoretical  task 
required  knowing  all  constituent  parameters  of  the  ZT  coefficient,  namely  the  Seebeck 
coefficient  S,  the  electrical  conductivity  a,  and  the  thermal  conductivity  k.  In  the  process  of 
developing  the  theory  of  the  Seebeck  coefficient  for  this  material,  a  new  transverse 
thermoelectric  phenomenon  at  zero  magnetic  field  was  discovered  which  we  call 
"anisotropic  ambipolar  transverse  thermoelectrics"  (AATT).  This  discovery  required  deeper 
theoretical  understanding  of  transverse  thermoelectric  phenomena  in  general,  with  the 
accompanying  realization  that  now  6  transport  coefficients  would  be  required  to 
characterize  this  material:  Sxx,  Syy,  kxx,  Kyy,  and  pxx,  pyy,  where  the  resistivity  tensor  p  is  the 
inverse  of  the  conductivity  tensor.  It  was  determined  that  S  and  p  tensors  could  be 
calculated  theoretically,  butK  would  have  to  be  measured  experimentally,  and  that  the 
single  component  Kyy  would  be  sufficient  to  estimate  device  performance.  In  the  end,  the 
AATT  behavior  of  the  T2SL  was  determined  to  give  a  transverse  ZiT  =  0.026  at  room 
temperature  which  is  able  to  produce  8  K  of  thermal  cooling  in  exponentially  tapered  10 
micron  thick  films,  and  can  yield  a  new  microscale  geometry  for  infrared  pixel  cooling.  It 
was  also  determined  that  wide  superlattice  thicknesses  which  have  a  topological-insulator¬ 
like  crossing  of  conduction  and  valence  bands  should  provide  greater  ZiT  values  at 
cryogenic  temperatures. 
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FIGURE  of  TYPE  II  SUPERLATTICE:  (left)  Band-diagram  llustration  of  type  II  superlattice  of 
InAs  (hosting  electrons)  and  GaSb  (hosting  holes).  Note  that  the  heavy-mass  holes  are 
confined  to  their  layer,  whereas  the  light-mass  electrons  can  freely  tunnel  between  layers. 
Proper  choice  of  chemical  potential  leads  to  a  semiconductor  with  net  n-type  conduction 
out-of-plane  and  net  p-type  conduction  in-plane,  (right)  Note  that  in  the  presence  of  a 
diagonal  electric  field  E  at  45  degrees,  the  net  charge  current  will  flow  parallel  to  the 
electric  field,  yet  the  net  particle  current  which  carries  heat  will  flow  orthogonal  to  the 
electric  field. 


The  remainder  of  this  Final  Report  is  therefore  broken  down  in  terms  of  five  principle  tasks: 

I.  Description  of  the  new  transport  equations  of  AATT  materials 

II.  Experimental  determination  ofKyy 

III.  Theoretical  determination  of  S  and  p  tensors  for  a  given  lattice  period 

IV.  Optimization  of  transverse  figure  of  merit  ZiT  over  all  lattice  periods  with  a 

description  of  a  new  device  geometries 

V.  Description  of  thin-metal-film  thermocouple  technology  needed  to  characterize 

thermoelectric  properties  of  T2SL  films  in  the  future 


3 


I.  Transport  equations  for  a  new  thermoelectric  phenomenon:  Ambipolar 
Anisotropic  Transverse  Thermoelectricity  (AATT) 

In  spite  of  the  prominent  success  of  standard  longitudinal  thermoelectrics,  there  remain  regimes 
of  operation  for  which  less  conventional  transverse  thermoelectrics  may  prove  advantageous. 
Longitudinal  thermoelectrics  are  extrinsically  doped  p-  or  n-type  with  heat  flow  parallel  or  anti¬ 
parallel  to  the  electrical  current,  respectively,  optimized  according  to  the  figure  of  merit  ZT  = 
Sxx^/pxxKxx  ,  with  most  recent  progress  above  room  temperature  using  nanoscale  phonon  scatters 
to  limit  the  thermal  conductivity  kxx  ,  while  maintaining  low  electrical  resistivity  pxx  and  large 
Seebeck  coefficient  Sxx.  Transverse  thermoelectrics,  on  the  other  hand,  have  a  large  off-diagonal 
component  Sxy  in  the  Seebeck  tensor  which  causes  heat  to  flow  orthogonal  to  the  electrical  current 
with  transverse  figure  of  merit  ZiT  =  Sxy2T/pxxKyy .  Whereas  the  use  of  longitudinal 
thermoelectrics  at  low  temperatures  is  restricted  since  the  dopants  freeze  out  at  cryogenic 
temperatures,  transverse  thermoelectric  phenomena  are  intrinsic  since  both  electrons  and  holes 
transport  heat  so  there  is  no  dopant  freeze-out.  And  whereas  the  minimum  temperature  for 
multistage  longitudinal  Peltier  coolers  is  limited  by  electrical  and  thermal  contact  resistances  at 
each  stage,  for  transverse  thermoelectrics  the  heat  and  electrical  current  flow  cross-sections  are 
orthogonal,  and  exponentially  tapered  devices  function  as  infinite-stage  Peltier  refrigerators, 
predicted  to  cool  to  arbitrarily  low  temperatures. [l]  Regarding  device  scaling,  the  fabrication  of 
microscale  longitudinal  thermoelectric  devices  can  be  quite  challenging.  [2]  On  the  other  hand, 
transverse  thermoelectrics  require  only  a  single  thermoelectric  layer,  so  devices  can  be  fabricated 
down  to  the  nanoscale.  Finally,  with  longitudinal  thermoelectrics  both  p-  and  n  -doped  materials 
need  to  be  separately  optimized,  but  with  transverse  thermoelectrics,  only  one  material  solution  is 
required. 

In  transverse  thermoelectric  phenomena  the  directional  symmetry  of  the  Seebeck  tensor  needs  to 
be  broken.  Thermomagnetic  effects  such  as  the  Ettingshausen  effect  require  an  external  magnetic 
field  to  break  the  time-reversal  symmetry  thereby  introducing  off-diagonal  terms  in  the  Seebeck 
tensor  and  generating  transverse  heat  flow.3  Composite  transverse  thermoelectrics  have  also  been 
proposed  which  have  a  structural  asymmetry  imposed  by  stacking  macroscopic  millimeter-sized 
slabs,4  alternating  a  semiconductor  with  large  Seebeck  coefficient  and  a  (semi)metal  with  large 
electrical  and  thermal  conductivity.  This  composite  can  induce  transverse  heat  flow  when  the 
current  is  skewed  with  respect  to  the  layers.  However,  the  usage  of  a  high  magnet  field  of  1.5  T  in 
the  former  case  limits  the  practical  application  of  the  Ettingshausen  effect,  and  the  stacking  of 
macroscopic  millimeter-scale  slabs  in  the  latter  case  does  not  permit  small  scale,  compact  devices. 

This  Final  Report  describes  the  band-engineering  strategy  which  we  call  anisotropic  ambipolar 
transverse  thermoelectrics  (AATT),  whereby  the  anisotropic  electron  and  hole  mass  tensors  of  a 
narrow  gap  superlattice  result  in  orthogonal  electron  and  hole  currents  and  a  large  transverse 
Seebeck  coefficient.  The  general  requisite  properties  for  AATT  are  outlined.  Transport  equations 
based  on  an  electron-hole  two-band  model  define  the  optimal  angle  of  electric  current  for 
inducing  the  maximum  transverse  figure  of  merit  ZlT.  The  concept  of  a  crossover  electric  field  Ec 
in  exponentially  tapered  devices  is  introduced  and  shown  to  distinguish  thin  and  thick  samples. 

The  equations  for  thermoelectric  transport  in  an  exponentially  tapered  cooler  are  also  explored, 
which  in  theory  can  yield  an  arbitrarily  low  cooling  temperature  with  sufficiently  thick  layers 
allowing  even  poor  ZiT  materials  to  induce  large  thermal  differences. 

We  begin  with  a  study  of  how  anisotropic  electron  and  hole  two-band  conduction  can  lead  to  the 
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desired  Seebeck  anisotropy:  an  ambipolar  Seebeck  tensor  which  is  n-type  in  one  direction  and  p- 
type  perpendicular.  Consider  the  general  case  of  a  semiconductor  near  intrinsic  doping  with 
anisotropic  electron  and  hole  conductivity  tensors  om  op  and  isotropic  Seebeck  tensors  s„,  sp 
where  x  and  y  axes  align  with  the  principle  axes  of  anisotropy: 
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where  s„<  0,  sp>  0.  The  total  conductivity  tensor  2  and  total  resistivity  tensor  P  are  related  as  2 
=  P  1  =  on  +  op.  The  total  Seebeck  tensor  is  the  weighted  sum  of  the  band  Seebeck  tensors, 

s  =  {<Tp  +  <Tn)-\apsp  +  <Tnsn)  .  (2) 


Whereas  single-band  Seebeck  tensors  s„,  are  typically  isotropic,  the  conductivity  tensors  om  op 
can  be  strongly  anisotropic,  and  Eq.  (2)  allows  one  to  exploit  this  anisotropic  conductivity  to 
weight  the  resulting  total  Seebeck  tensor  to  opposite  signs  of  Seebeck  effect  for  orthogonal 
directions.  Defining  a  small  parameter  §  as  the  ambipolar  conductivity  ratio  for  a  given  direction, 
we  shall  define  the  v-direction  as  dominated  by  p-type  conduction  ^x  =  on>xx=_op>xx  <  1 ,  and  the  y 
-direction  by  n-type  ^  =  op>yyl  on>yy<  1.  The  total  Seebeck  tensor 


Sp,XX  0 
0  *^n,yy 


(3) 


has  the  following  positive  and  negative  matrix  elements,  respectively,  for  small  %  <  sp  sn  /  (s» 
s2) 
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(4) 


To  describe  the  transverse  thermoelectric  behavior,  we  define  longitudinal  and  transverse 
transport  axes  a  and  b ,  respectively.  Let  the  unit  vector  a  =  cos^x  +  sin^v  define  the  direction  of 
current  flow  J  =  *^aa  at  an  angle  +9  relative  to  the  x  -axis,  with  the  transverse  unit  vector 
b  =  —  sin0x  —  cos0y?  as  shown  in  the  insert  of  Fig.  4.  The  total  Peltier  heat  flux  density  becomes 
Q„  =  nj  =  (TS)J  with  longitudinal  and  transverse  components, 

Qn,a  =  Qn  a=  (Sp,xxcos20  +  Sn,yysm20)T Ja  (5) 

Qn,b  =  Qn  •  b  =  —  (SPyXX  —  Sn^yy)  cos0sin0 TJa.  (6) 

Note  from  Eq.  (6)  that  the  transverse  Peltier  effect  can  be  maximized  by  independently 
maximizing  the  magnitude  of  Sp>xx  and  Slhyy  . 

The  total  heat  flux  density  Q  =  Qn  -  kcVT  will  include  contributions  from  both  Peltier  and 
thermal  conduction  effects,  where  Kc  with  the  explicit  superscript  c  defines  the  open-circuit 
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7  _  r  -  X7  T  =  — h 

thermal  conductivity  tensor  at  J  =0.  Assuming  J  —  and  db  ,  Kooi  et  al.  have  shown 

that  the  longitudinal  electric  field  component  is  constant  everywhere.  Assuming  the  thermal 
gradient  VT  is  everywhere  perpendicular  to  the  heat  sink  in  the  b  =  0  plane,  the  heat  flux 
component  Q/,  depends  only  on  b. 

The  longitudinal  current  and  transverse  heat  flow  are 

T  —  -J—fP  Egb  dT 

a  *-■  a  Paadb 

jrr 

Qt=  Tj£Ea  -(i+Z^kU-jg, 

ZbgT  =  ^ al>  ^ba  ^ 

where  the  transverse  figure  of  merit  is  a  PaaKw>  . 

Under  steady  state,  the  current  and  energy  flux  densities  have  zero  divergence  V*  J  =  0  and  V#(Q 
+  jiJ  )  =  0,  where  p  is  the  electrochemical  potential,  and  -Vp  =  E  defines  the  electric  field.  The 
intuitive  result  is  that  Joule  heating  EJa  from  the  longitudinal  current  density  sources  a 
transverse  heat  flux  density  Qa  : 


(7) 
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Equations  (7)-(9)  define  the  differential  equation: 
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Assuming  a  limited  temperature  range  with  constant  thermoelectric  coefficients,  we  get: 

(ii) 


P  _  ,  Ea  dT .2  1  +  ZbgT  (PT 

Zba  dV 


Notice  this  equation  differs  from  that  for  the  magnetic  field-induced  Ettingshausen  effect,  since  in 
that  case  Sab  and  Sba  have  opposite  signs  and  all  dT/db  terms  cancel;  whereas  in  AATT  the  Sab  and 
Sba  are  equal,  preserving  the  dT/db  term  in  Eq.  (1 1). 

The  optimal  angle  0±  which  maximizes  ZbaT  (0)  defines  the  optimal  ZiT 

COS20±=  - !  1  (12) 


1i/  K3Lil(K* 

V  pvvfp** 

Z,T  =  ZM)T  =  {JZZ:Z%Z\-  (13) 


Note  that  the  angle  0±  is  independent  of  the  Seebeck  anisotropy,  and  approaches  0X  =  jt/4  as  the 
thermal  conductance  anisotropy  matches  the  resistance  anisotropy  Kyy/Kxx  =  pyy/pxx  .  To  optimize 
the  electronic  band  structure  it  is  useful  to  define  a  power  factor  PF±±  from  Eq.  (13)  under  the 
assumption  of  isotropic  k  where  S  and  p  tensors  can  be  readily  calculated  with  simple  scattering 
assumptions: 
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II.  Experimental  determination  of  thermal  conductivity,  Kyy 

The  two  different  T2SLs  studied  in  this  work  were  grown  by  molecular  beam  epitaxy  on  GaSb 
substrates.  Because  this  material  was  designed  for  p-i-n  detectors,  it  consists  of  a  0.5  pm  GaSb  p+ 
(~  1018  cm3)  buffer  layer,  followed  by  a  0.5  pm  T2SL  p+  (~  1018  cm3 )  region,  a  2  pm  undoped 
T2SL  layer,  a  0.5  pm  T2SL  n+  (~  1018  cm3)  region,  and  a  10  nm  Si  doped  InAs  n+  capping 
layer.  T2SL-1  is  composed  of  12  monolayers  (ML)  of  InAs  and  8  ML  of  GaSb  per  period.  T2SL- 
2  is  composed  of  19  ML  of  InAs  and  18  ML  of  GaSb  per  period. 

Because  of  the  similarities  in  composition,  the  thermal  conductivity  in  this  T2SL  should  be 
representative  of  the  thermal  conductivity  in  generic  T2SL  device  layers,  including  the  T2SL 
cladding  layers  of  InAs/AlSb  and  T2SL  active  layers  of  InAs/GalnSb. 

Following  standard  3co  sample  preparation,  the  thermal  conductivity  of  a  sample  including  the 
layer  of  interest  is  measured  relative  to  a  reference  sample  without  the  layer.  Thus  the  T2SL  is 
wet-etched  away  from  half  the  sample  with  a  solution  of  citric  acid  and  phosphoric 
acid  plus  peroxide. 


Fig.  1  Temperature  change  AT  vs.  frequency  /  measured 
with  the  3a;  method  for  the  T2SL  +  SiC>2  +  GaSb  substrate 
(filled  symbols)  and  for  SiC>2  +  GaSb  substrate  (open  sym¬ 
bols),  (a)  at  15  K  and  300  K,  and  (b)  at  80  K  and  212  K. 
Inset:  picture  of  the  sample,  with  left  half  etched  to  sub¬ 
strate,  and  right  half  unetched.  Two  Au  heater-thermometer 
lines  are  deposited  on  each  side,  with  current  contact  pads  at 
top/bottom  and  voltage  contacts  in  between. 


An  insulating  Si02  layer  of  150  nm  is 
deposited  using  plasma  enhanced 
chemical  vapor  deposition  (PECVD)  to 
prevent  an  electrical  short  circuit 
through  the  conducting  substrate.  Then 
200  nm  thick  gold  heater-thermometer 
filaments  are  deposited  atop  a  3  nm  Ti 
adhesion  layer  using  e-beam 
evaporation  on  both  etched  and 
unetched  regions.  The  filament  shown 
in  the  inset  of  Fig.  1(b)  is  3.6  mm  long 
and  30  pm  wide,  much  wider  than  the  3 
pm  T2SL  thickness  so  that  the  heat 
flow  through  the  T2SL  obeys  the  one¬ 
dimensional  thermal  diffusion  equation. 


According  to  the  3oo  method, [5,6]  by 

sending  a  current  excitation  at  frequency  co  through  the  gold  filament,  heating  power  P  is  induced 
at  frequency  2oo,  leading  to  thermal  variation  AT  at  the  same  frequency  and  ultimately  a  4-point 
voltage  V3(0  at  frequency  3co  proportional  to  the  temperature  difference  via  the  thermal  coefficient 
of  Au.  The  thermal  difference  A T=  ATsub  (co  )  +  A Tf  can  have  a  frequency  dependent  substrate 
contribution  and  a  frequency  independent  thin-film  contribution, 


4Tsub(u,)=  — Liln(2a;)  +  C 

»7TKgub  ^ 


at,  =  -2- 

IW 
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where  Ksub  is  the  thermal  conductivity  of 
the  GaSb  substrate,  w  and  /  the  width 
and  length  of  the  Au  filament  heater, 
and  C  is  a  constant  offset  which 
includes  the  Si02  layer  contribution.  For 
the  T2SL  film,  Kf  and  t  are  the  film 
thermal  conductivity  and  thickness, 
respectively. 

We  measured  the  sample  in  an  Oxford 
variable  temperature  insert  (VTI) 
helium  gas  flow  cryostat  from  300  K 
down  to  13  K,  using  standard  lock-in 
techniques.  Since  there  are  background 
3 co  voltages  from  the  lock-in  power 
source  and  nonlinear  components  in  the 
measurement  circuit  and  in  the  lock-in 
A-B  input  channels,  we  also  measured  a 
reference  background  3co  signal  with  a 
low  thermal  coefficient  resistor  of  equal 
resistance  as  the  gold  filament. 

By  measuring  the  slope  of  the  temperature  difference  AT  as  a  function  of  log-frequency,  we  can 
deduce  the  substrate  thermal  conductivity  and  thereby  confirm  the  reliability  of  our  measurements. 
In  Fig.  2,  we  compare  our  measured  GaSb  substrate  thermal  conductivity 
(solid  circles)  with  previously  published  GaSb  bulk  thermal  conductivity  (open  circles)[7,8], 
indicating  excellent  agreement. 

The  T2SL  cross-plane  thermal  conductivity  Kf  is  deduced  from  the  thermal  difference  A  7> , 
plotted  with  solid  triangles  in  Fig.  2.  Previous  studies  of  superlattice  thermal  conductivities  have 
shown  that  superlattice  thermal  conductivity  is  significantly  reduced  from  those  of  their 
constituent  bulk  materials,  resulting  from  phonon  interface  scattering,  reduced  group  velocity 
from  modified  phonon-dispersion  relation,  and  strain-relaxation  induced  high  density  of 
dislocations  [9,10,n].  The  T2SL  value  is  reduced  by  2  orders  of  magnitude  compared  with  the  bulk 
substrate  thermal  conductivities  for  GaSb  bulk  (open  circles).  We  note  that  the  suppression  is 
much  greater,  up  to  as  much  as  3  Vi  orders  of  magnitude,  when  compared  with  InAs  bulk  (open 
squares). 

III.  Theoretical  determination  of  S  and  p  tensors 

Calculations  of  Seebeck  S  and  resistivity  p  tensors  for  the  specific  case  of  type  II  broken-gap 
InAs/GaSb  superlattices  (T2SL)  are  shown  to  give  promising  values  at  room  temperature  for  use 
as  nanoscale  transverse  thermoelectric  refrigerators  and  generators.  Wide-period  T2SL’s  are 
proposed  to  be  likely  candidates  for  low-temperature  cryogenic  thermoelectrics. 

We  can  apply  the  power  factor  optimization  of  Eq.  (14)  to  the  InAs/GaSb  type  II  broken-gap 
superlattice  (T2SL),  and  identify  it  as  a  likely  candidate  AATT  thermoelectric  with  its  tunable 
anisotropic  electrical  conductivity  tensor  and  tunable  band  gap.  The  T2SL  is  a  well-studied 
material  for  low-dark  current  infrared  detectors  and  emitters.  Because  of  the  broken-gap 
alignment  at  the  heterointerface,  the  valence  band  edge  of  GaSb  lies  above  the  conduction  band 
edge  for  InAs.  Thus  by  creating  superlattices  of  varying  layer  width,  the  quantum  confinement 


Fig.  2  Cross-plane  thermal  conductivity  for  T2SL  samples 
(solid  triangles).  Measured  GaSb  substrate  thermal  conduc¬ 
tivity  is  shown  in  solid  circles.  Published  data  for  bulk  ther¬ 
mal  conductivity  of  InAs  [19]  and  GaSb  [20]  are  shown  for 
comparison  in  open  squares  and  open  circles  respectively.  The 
solid  line  is  a  fitting  polynomial  for  the  average  value  of  T2SL 
thermal  conductivity  described  in  the  text. 
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energy  can  produce  a  spatially  indirect  semiconductor  with  arbitrarily  small  energy  gap. 
Furthermore,  since  the  heavy  hole  mass  prohibits  the  holes  in  GaSb  from  tunneling  through  the 
InAs  layers,  no  hole  current  flows  perpendicular  to  the  superlattice,  and  the  T2SL  hole 
conductivity  is  maximally  anisotropic.  One  can  therefore  adopt  a  3D  ellipsoidal  effective  mass 
approximation  for  the  electrons  and  a  2D  circular  effective  mass  approximation  for  holes  to 

( a  z  t + Ei  ^ 

V  *Bt  ) 

(*©*) 

CTp’xx  =  ^2(fcBT)fl+2F‘+*  (a^t) 

^p*yy  —  0  (15) 

Sn  =  — [s  +  5/2  +  (Eg  -  fi)/kgT] 
e 

Sp  —  - [s  +  5/2  +  fi/kBT] 

e 


determine  the  conductivity  and  Seebeck  tensor  components: 

=  ^^rV^(kBTy+5/2FV2+s 

y^(A-BT)a+5/2F3/2+s 


where  mx  represents  the  effective  mass  component  of  the  band  of  interest  n  or  p  in  the  direction 
of  interest  x  or  y,  ji  is  the  chemical  potential  measured  relative  to  the  valence  band-edge,  and  F  is 

n  /  r cSD  mj-  | 

the  Fermi-Dirac  integral  y'r^/  =  ./o  ^  i+e*  .  The  scattering  time  x  =  yEs  where  E  is  the  kinetic 
energy  of  the  carrier.  Since  the  scattering  is  dominated  by  interface  scattering,  we  assume  the 
power  factor  5  =  0.  Then  y  can  be  calculated  from  equation  y  =  \ic  m*/e  ,  where  jic  is  the  carrier 
mobility,  m*  is  the  carrier  effective  mass  and  e  is  the  electron  charge. 


InAs:  7.2701  nm 
GaSb:  7.0103  nm 


InAs:  7.8759  nm  InAs:  9.0876  nm 

GaSb:  7.0103  nm  GaSb:  8.5342  nm 


k  (1/Angstrom) 


k(l/Angstrom) 


k(l/Angstrom) 


FIG  3.  Band  simulations  of  the  dispersion  relation  of  a  superlattice,  show  semiconducting 
behavior  (left),  zero  gap  (center),  and  inverted  semimetallic  behavior  (right).  Note  the  small 
anticrossing  gaps  on  the  right  form  a  highly  degenerate  a  ring-shape  in  momentum  space,  which 
functions  as  a  large  valley-degeneracy,  good  for  enhancing  the  Seebeck  coefficient  at  cryogenic 
temperatures. 

These  T2SL  materials  should  have  an  enhanced  low-temperature  performance  in  wide -period 
superlattices.  Simulations  have  been  performed  on  a  superlattice  structure  in  Fig.  3,  showing  that 
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by  varying  the  well  thickness  it  is  possible  to  transition  the  k  =  0  band  gap  from  semiconducting, 
to  zero  gap,  to  semimetallic,  just  as  with  the  double-well  case.  Note  that  the  small  anticrossing 
gap  in  the  “semimetallic”  case  will  actually  result  in  a  large  density  of  states  analogous  to  a 
valley-degeneracy,  and  a  small  gap  —  both  hallmarks  of  a  good  low-temperature  thermoelectric 
material.  We  propose  that  this  discovery  identifies  an  interesting  avenue  of  future  research  for 
cryogenic  transverse  thermoelectrics. 

IV.  Optimization  of  transverse  figure  of  merit  ZiJ  and  new  device  geometries 

We  calculate  the  band  structure  of  the  T2SL  with  the  nextnano3  full-band  envelope  function 
method.  The  figure  below  shows  on  the  left  how  the  phase  space  of  superlattice  layer  thickness, 
^GaAs  vs.  dinAs  is  mapped  out.  Each  point  in  this  phase  space  corresponds  to  a  band  structure 
calculation,  as  shown  on  the  right.  Going  from  narrow-  to  wide-periods,  one  encounters  a  tunable 
gap  which  goes  from  semiconducting  gap  to  zero  gap  to  semi-metallic.  All  transport  parameters 
derived  the  band  structure,  namely  the  energy  gap  and  effective  masses,  are  all  calculated  for 
every  point  on  the  dG aAs  vs.  di^  phase  space,  and  the  optimal  power  factor  is  determined. 


Type  II  Superlattice 


E 

.o 

09 

■ 

O 


£ 

*o 

I 


NextNANO 3:  Full-band 
envelope-function 
simulation  InAs/GaSb 
superlattice  dispersion 


Positive  gap: 
Semiconductor 


Zero  gap 


Semimetal 


Applying  the  results  to  Eqs.  (14)  and  Eq.  (15),  one  can  maximize  the  power  factor  PF±  at  300  K. 
The  result  yields  optimal  with  layer  thicknesses  dh ^  =  7.88  nm  and  dG asb  =  3:  96  nm.  The 
electron  and  hole  effective  masses  are  mn>x  =  0:  0303  m0  ,  mn>y  =  0.0268  mo  and 
mp>x  =  0:  0474  mo  where  mo  is  the  free  electron  mass.  The  energy  gap  is  Eg  =  28.3  meV,  and  the 
maximum  power  factor  occurs  at  chemical  potential  q  =  2.43  meV.  Under  the  assumption  that  the 
thermal  conductivity  is  isotropic,  the  optimal  angle  is  0±  =  32°,  which  means  that  most  of  the  heat 
flow  will  be  perpendicular  to  the  superlattice,  so  the  thermal  conductivity  is  dominated  by  the 
out-of-plane  component  experimentally  measured  in  the  previous  section  to  be  k  =  4  W/m*K.  The 
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optimal  ZT  is  therefore  estimated  to  be 
Z±T  =  0.026  at  room  temperature. 

The  maximum  AT  is  reached  when  the 
heat  load  at  the  cold  side  is  zero.  So 
the  boundary  conditions  for  Eq.  (11) 
are  Th  =  300  K  and  Qm  (Tc  )  =  0.  The 
temperature  profile  is  shown  in  Fig.  4 
for  a  cross  section  thickness  of  b  =10 
qm  of  a  superlattice,  yielding  a  total 
temperature  difference  ATmax  =  4  K. 
We  also  solve  the  temperature  profile 
with  ZTw  =0.1  or  ZTw  =  0.5  with  the 
maximum  temperature  differences  of 
14  K  and  56  K,  respectively,  and  we 
compare  with  both  longitudinal  and 
Ettingshausen  cooling  for  the  same  ZT 
values  as  shown  in  Fig.  4.  AATT 
transverse  thermoelectric  cooling 
(solid  line)  is  slightly  better  than  the 
longitudinal  thermoelectric  cooling 
(dotted  line),  and  the  thermal  profile  is 
distinctly  different  from  the  Ettingshausen  effect  for  a  hypothetical  material  with  the  same  ZT 
value  (dashed  line). 

More  dramatic  improvements  in  cooling  power  can  be  expected  for  transverse  thermoelectrics  by 
exploiting  tapered  geometries.  It  has  long  been  known  that  for  longitudinal  Peltier  coolers,  a 
cascade  of  ideal  lossless  thermoelectric  junctions  with  arbitrarily  small  ZT  and  multiple 
stages  should  cool  to  arbitrarily  low  temperatures.  The  problem  with  longitudinal  thermoelectrics 
is  that  the  electrical  and  thermal  contact  resistance  at  the  interfaces  between  stages  limits  the 
practical  implementation  to  at  most  6  stages12,13  ,  and  the  maximum  AT  is  limited 
to  ViZTq2  for  a  single  stage.  In 
transverse  thermoelectrics,  since 
the  current  and  heat  flow  are 
perpendicular  to  each  other, 
infinite  cascading  can  be 
effectively  implemented  by 
making  the  device 
exponentially  tapered  as  it 
approaches  the  cold  load  side, 
as  in  the  insert  of 
Fig.  5(b).  Following  the 
analysis  of  Kooi  et  al.  for  the 
Ettingshausen  effect,  this 
exponential  tapering  of  the 
device  in  the  perpendicular  c- 
direction  adds  the  term  of  Qb 
d(ln  c)/Ab  to  Eq.  (9).  Assuming 
exponentially  tapering  c  =  c0e 
,  where  L  is  the  length  scale 
for  exponential 
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FIG.  5:  Temperature  profile  solved  with  boundary  conditions 
Th(0)  =  300  K,  Qin(Tc)  =  0  at  different  values  of  electric 
field  gives  the  ATmax  at  different  device  thicknesses.  The 
coordinates  of  the  squares  and  circles  identify  the  minimum 
Tc  for  that  layer  thickness,  (a)  For  the  rectangular  shape 
device,  the  minimum  cooling  temperature  is  independent  of 
the  layer  thickness,  but  weaker  electric  fields  require  thicker 
layers.  The  electric  fields  Ea  shown  are  at  step  intervals  of 
0.005  V /fin l.  (b)  For  the  exponentially  tapered  device,  ATmax 
can  be  arbitrary  large  with  increasing  layer  thickness. 


FIG.  4:  Temperature  profile  for  Ettingshausen  (dashed  line), 
transverse  thermoelectric  (solid  line)  and  longitudinal  ther¬ 
moelectric  (dotted  line)  cooling  for  various  ZT  values.  Heat 
sink  is  at  6  =  0,  Th  =  300  K,  and  Qin  ( b  =  10/xm)  =  0.  Elec¬ 
tric  field  is  chosen  to  maximize  AT.  Z  is  assumed  temperature 
independent,  and  Sab  =  -400  /xV/K. 


11 


tapering,  Eq.  (1 1)  for  AATT  materials  becomes: 


Fig.  5  shows  the  temperature  profile  solved  with  different  device  thicknesses  at  corresponding 
values  of  electric  fields  for  both  a  rectangular  cooler  without  tapering,  and  an  exponentially 
tapered  device.  Note  that  for  the  rectangular  cooler  in  Fig.  5(a),  the  ATmax  are  independent  of 
layer  thickness,  whereas  in  Fig.  5(b)  the  exponential  tapering  gives  an  ever  ATmax  with 
increasing  sample  thickness.  Assuming  the  same  exponential  tapering  factor  as  was  demonstrated 
in  bismuth  semimetals14  b  =  L  =  2.77,  the  AATT  cooling  becomes  ATmax  =  8  K.  This  is 
competitive  with  recent  experimental  results  in  on-chip  cooling  by  significantly  more  complex 
device  structures  made  of  superlattice-based  longitudinal  thin-film  thermoelectrics.  Solving  Eq. 
(16)  for  extrema,  we  find  a  critical  crossover  field  Ec  =-  SahTH/L  for  exponentially  tapered 
transverse  thermoelectric  devices.  At  Ea  =  Ec ,  the  AATT  cooling  is  balanced  with  Joule  heating 
at  each  cross-section  perpendicular  to  b,  and  the  temperature  is  uniform  in  the  whole  device.  Ea  > 
Ec  is  required  for  a  local  maximum  in  AT7  to  exist,  wherefore  a  limit  exists  for  ATmax.  When  E a  < 
Ec ,  on  the  other  hand,  T{b )  decreases  monotonically  and  there  is  no  limit  for  AT(Z>),  which  means 
with  a  small  electric  field  and  a  sufficiently  thick  slab  of  material,  the  device  can  be  cooled  to 
arbitrarily  low  temperature.  In  reality,  the  material  transport  coefficients  are  temperature 
dependent,  which  limits  ATmax  .  As  shown  in  Fig.  5(b),  for  thin  devices  where  b=L  is  small, 

ATmax  are  given  by  Ea  >EC  shown  with  the  solid  lines.  Whereas  for  thick  devices  where  b=L  is 
large,  ATmax  are  given  by  Ea  <EC  shown  with  the  dashed  lines. 


|  Heat  current 


|  Electrical  current 


1+ 


FIGURE  of  Unconventional  transverse  thermoelectric  cooling  geometry.  An  anular  cathode  is 
placed  atop  a  T2SL  transverse  thermoelectric,  fed  by  a  central  anode  below,  surrounding  an 
infrared  pixel  which  is  to  be  cooled.  As  electrical  current  (black  arrows)  flows  between  the 
electrodes,  the  heat  current  flows  opposite  the  vertical  current,  but  parallel  to  the  lateral  current. 
In  between,  the  transverse  thermoelectric  effect  causes  currents  at  45  degrees  to  pull  heat  away 
from  the  central  area.  At  every  point,  the  pixel  is  being  cooled  by  the  electrical  current. 
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IV.  Thin-metal-film  thermocouple  technology 


(a) 


(b) 


Fig.  6:  (a)  Conceptual  drawing  of  bimetallic  thermocouple  layers,  where 
metal  A  represents  the  first  metal  deposited,  while  metal  B  represents  the 
second.  Heat  was  applied  on  backside  of  substrate  behind  the  hot 
thermocouple  at  7h.  whereas  the  reference  cold  thermocouple  was  at  7c.  In  the 
real  device  the  left  and  right  leads  of  Metal  A  wrap  around  to  form  voltage 
leads.  A  resistance  temperature  dependent  (RTD)  thermistor  is  epoxied  on  top 
of  each  bimetallic  junction  calibrating  the  absolute  thermocouple 
temperatures,  (b)  Microscope  photograph  of  50  nm  thin  film  thermocouple  of 
Cr  (Metal  A)  and  Ni  (Metal  B).  on  a  semi-insulating  GaAs  substrate.  Seebeck 
voltage  is  measured  between  leads  at  Tq. 


Thin  film  thermocouples  are  tested  for 
eventual  deposition  on  Type  II 
superlattices.  The  thermocouples  can 
be  placed  on  the  superlattice  with  an 
intervening  oxide  layer  for  electrical 
isolation.  These  thermocouples  will 
calibrate  thermal  gradients  across  or 
along  the  superlattice  so  that  the 
Seebeck  coefficient  and  thermal 
conductivity  of  the  superlattice  itself 
can  be  calibrated.  The  best  choice  for 
metals  was  determined  to  be  nickel 
and  chromium,  producing  a  41  mV/K 
thermocouple  signal.  Note  that  once 
transverse  thermoelectric  behavior  had 
been  identified,  the  efforts  towards 
Seebeck  and  conductivity  calibrations 
were  diverted  to  the  thermal 
conductivity  measurements  described 
in  Section  II,  above. 


The  bimetallic  thin-film  thermocouples  were  deposited  on  blank  semiconductor  substrates  to  test 
their  functionality.  This  newly  developed  method  for  measuring  temperature  gradients  can  be 
readily  implemented  with  existing  fabrication  processes.  We  first  identify  specific  thin-film 
metals  with  a  sufficiently  large  Seebeck  coefficient  and  then  demonstrate  reproducible  process 
technology  on  semiconductor  substrates.  Other  impacts  of  this  work  include  bimetallic 
thermocouple  sensors  for  integration  into  microchip  fabrication  to  achieve  thermal  accuracies 
within  a  fraction  of  a  degree  Celsius.  Samples  produced  in  this  research  have  demonstrated 
reproducible  scaling  factors  for  both  inter-  and  intra-batch,  as  required  for  industry  applications. 
Another  possible  advantage  of  this  technique  is  that  it  is  entirely  passive,  allowing  in  principle  for 
thermal  mapping  of  test  samples  and  microprocessor  chips  with  arbitrary  resolution  and  zero 
power  overhead ,  allowing  cryogenic  calibrations  with  no  load  heating. 


The  two  metals  tested  here  in  a  thin  film  thermocouple  are  chromium  Cr  and  nickel  Ni.  The 
thermocouple  discussed  here  is  shown  in  Fig.  6(a),  with  metal  A  representing  the  first  metal 
deposited  Cr,  and  metal  B  representing  the  second  metal  deposited  Ni. 

The  Seebeck  effect  arises  when  a  temperature  difference,  TH  -  Tc,  occurs  between  two  junctions 
of  dissimilar  conductors.  An  open-circuit  voltage  proportional  to  the  junction  temperature 
difference  is  thereby  induced  at  the  leads  independent  of  the  lead  temperature  To, 


Vhc  ~  Tc) 


(17) 


Here  Vhc  represents  the  Seebeck  voltage  between  the  leads.  The  temperatures  Tu  and  Tc  are  the 
temperatures  located  at  the  hot  and  cold  thermocouple  junctions,  respectively,  as  noted  in  Fig. 
6(b).  The  temperatures  are  calibrated  with  small  area  (2.0  mm  x  3.0  mm)  resistance  temperature 
dependent  (RTD)  thermistors  which  are  epoxied  with  heat-conducting  silver  epoxy  on  top  of  the 
metal  junctions  as  shown  in  Figure  6(a).  a  is  the  Seebeck  coefficient  of  the  thermocouple.  The 
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bulk  Seebeck  coefficient  aB  expected  from  literature  can  be  determined  by  taking  the  difference 
of  Seebeck  coefficients  of  two  dissimilar  conductors.  Cr  and  Ni  were  chosen  for  their  relatively 
high  bulk  Seebeck  coefficients  of  +21.8  and  -19.5  jiV/°C  respectively,  to  produce  a  predicted  aB 
=  41.3  |iV/°C.  For  this  study,  pure  metals  were  chosen.  Though  the  various  alloys  of  Ni/Cr  hold 
promise  for  thermocouple  applications,  alloys  are  also  known  to  fractionate  upon  evaporation, 
possibly  making  Seebeck  calibration  less  reliable  and  a  topic  for  future  study. 

Thermocouple  devices,  as  seen  in  Fig.  6(b),  were  fabricated  on  semi-insulating  (SI)  GaAs 
substrate.  The  GaAs  substrate  had  a  thickness  of  450  pm  and  a  high  bulk  resistivity  of  3.3  x  106 
Cl- m.  GaAs  was  chosen  as  the  semiconductor  substrate  both  for  its  high  resistivity  and  easy 
cleaving  into  narrow  strips.  Both  metal  layers  were  50  nm  thick,  made  by  using  positive 
photoresist  photolithography  and  electron-beam  evaporation.  To  remove  any  oxide  between  the 
two  metal  layers,  an  HC1  acid  dip  was  performed  immediately  before  the  second  metal 
evaporation.  The  long,  narrow  thermocouple  design  in  Fig.  6(b)  is  chosen  to  provide  large  linear 
thermal  gradients  for  a  given  heater  power,  allowing  for  better  calibration. 


Experiments  were  conducted  by  mechanically  clamping  the  cold  end  of  the  semiconductor 
substrate  in  between  two  blocks  of  aluminum,  which  functioned  as  a  cold  sink.  Gold  wires 
connected  the  leads  of  the  device  to  lines  of  indium  on  a  glass  slide.  Copper  wires  were  then 
connected  from  the  lines  of  indium  on  the  glass  slide  to  digital  multimeters.  To  simulate  an 
integrated  circuit  hot  spot,  a  soldering  iron  was  used  as  the  heater  element,  applied  to  the  backside 
of  the  GaAs  substrate  under  the  thermocouple.  The  heater,  set  at  room  temperature,  was  clamped 
into  place  to  make  contact  with  the  back  side  of  the  device  at  the  location  of  TH  as  seen  in  Fig. 
6(b).  Then  the  heater  power  was  gradually  increased  to  a  desired  temperature,  and  then  gradually 
cooled  back  down  to  room  temperature.  Hysteresis  artifacts  due  to  uneven  heat  flow  can  thus  be 
identified  and  avoided.  This  slow  heating  process  limited  the  potential  adverse  effects  of  rapidly 
heating  the  thermocouples  which  may  result  in  non-linear  thermal  gradients. 


Experimental  data  was  collected  for  AEHc,  TH,  and  Tc  for  both  increasing  and  decreasing 
temperature.  Because  Tc  itself  varies  with  the  heat  load,  any  non-linearities  in  the  Seebeck 

response  would  induce  false  hysteresis 
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in  the  thermally  cycled  raw  data.  Thus 
the  following  relations  were  used  to 
correct  for  small  thermal  variations  of 
Tc- 


-  (^H  TQ) 
Arco  =  (Tc  —  T0) 


dV 

AFH0=AFHC+—-A  tco 
dT 


(18) 

(19) 

(20) 


Fig.  7:  Parametric  plot  of  thin-film  Ni/Cr  thermocouple  Seebeck  response 
according  to  Eq.  (4).  with  absolute  temperature  at  the  top  axis  and  relative 
temperature  on  the  bottom  axis.  Traces  are  shifted  by  100  pV  for  clarity. 
Thermocouples  were  taken  from  different  batches  on  semi-insulating  (SI) 
GaAs  substrates  and  their  Seebeck  coefficients  were  read  by  the  slopes  on  the 
voltage  vs.  temperature  plot  above.  The  measured  Seebeck  coefficients 
reported  in  Table  1  was  smaller  the  bulk  value  by  the  factor  T  indicated. 


The  term  AfHo  represents  the  Seebeck 
voltage  as  a  function  of  ATHo  while 
AEhc  is  the  Seebeck  voltage  as  a 
function  of  A  The  determined  by  (17). 
The  term  dV/dT  is  the  linear  Seebeck 
coefficient  of  the  sample  for  small 
thermal  differences  near  ambient 
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temperature,  which  is  calculated  by  taking  the  average  aF  over  the  range  A  The  =  OK  to  A  The  = 

20  °C.  By  implementing  these  equations,  an  accurate  estimate  of  the  Seebeck  response  over  the 
full  temperature  range  can  be  made,  without  presupposing  linearity.  A7ho  and  AfHo  of  the  devices’ 
responses  are  plotted  in  the  parametric  plots  shown  in  Fig.  2. 


The  parametric  plots  of  the  Seebeck  response  were  observed  to  be  perfectly  linear  within 
measurement  noise,  so  that  a  single  Seebeck  coefficient  aF  of  the  thin  film  could  be  determined 
from  the  slope. 

aF=^ 

(21) 

Several  batches  of  50nm  GaAs  devices  were  fabricated  to  determine  the  design’s  reproducibility 
between  batches.  With  A  Tm  representing  the  absolute  temperature  change  from  room 
temperature  and  A  Vm  representing  the  corrected  Seebeck  voltage,  the  parametric  plots  in  Fig.  2 
exhibit  the  clear  reproducibility  of  the  Seebeck  response.  The  thermocouple  range  is  tested  up  to 
ATHo  =  60  °C  equivalent  to  an  absolute  temperature  of  TH  =  90  °C. 


Samples  produced  thin  film  Seebeck  coefficients  aF  smaller  than  the  bulk  literature  value,  so  a 
scaling  factor  T  was  defined, 


r  = 


aF  (22) 

The  statistics  of  the  devices  depicted  in  Fig.  7  are  presented  in  Table  I.  The  scaling  factors  for  the 
GaAs  #1,  GaAs  #2  and  GaAs  #3  samples  are  quite  similar. 


The  above  table  describes  the  important  result  of 
this  work,  namely  the  reproducible  factor  of  T  = 

3.98  scaling  between  the  50  nm  thin  film  Seebeck 
coefficient  and  the  bulk  Seebeck  coefficient.  The 
cause  of  such  thin-film  reduction  of  the  Seebeck 
effect  has  been  reported  experimentally  and 
investigated  theoretically,  and  is  shown  to  arise 
when  the  film  thickness  is  of  order  or  smaller  than 
the  electron  mean-free-path,  such  that  boundary  scattering  starts  to  compete  in  magnitude  with 
bulk  scattering,  typically  for  films  thinner  than  several  hundred  nanometers,  depending  on  the 
metal.  The  difference  between  the  thin  film  Seebeck  coefficient  and  that  of  bulk  material  is 
inversely  proportional  to  the  metal  film  thickness.  Another  possible  explanation  is  that  the  films 
may  incorporate  oxygen  during  deposition  resulting  in  a  reduced  Seebeck  coefficient,  though  this 
concern  is  more  relevant  for  rf  sputtered  films.  Between  the  two  consequent  electron-beam 
evaporations,  processing  may  result  in  a  thin  oxide  or  other  process  contamination  layer  which  in 
principle  could  also  reduce  the  Seebeck  coefficient.  The  likelihood  of  an  oxide  layer  is  reduced 
since  the  sample  is  dipped  in  HC1  acid  immediately  before  the  second  metal  evaporation. 
However,  even  if  the  scaling  factor  does  result  from  some  sort  of  extrinsic  layer,  these  conditions 
have  been  demonstrated  here  to  be  reproducible  in  several  attempts.  By  increasing  the  thickness 
of  each  of  the  metal  layers  the  Seebeck  coefficients  could  potentially  increase  closer  to  the  bulk 
literature  value.  At  higher  temperatures  up  to  700  °C,  Moore  et  al.  determined  that  nonlinearities 
may  enter  the  Seebeck  response  of  some  bulk  metals,  though  the  response  for  the  temperatures 
presented  here  is  demonstrably  linear.  Thus  for  a  given  layer  thickness,  we  have  shown  that  the 
Seebeck  coefficient  can  be  reproducibly  calibrated,  as  we  have  done  here  for  the  50  nm  case.  The 
thermocouples  can  be  scaled  down  and  routed  along  arbitrarily  shaped  paths  in  dedicated  metal 
layers  on  top  of  a  sample. 


TABLE  I 

Reproducibility  Statistics 


Substrate 

aF(Mvrc) 

r 

GaAs  #1 

10.19 

4.05 

GaAs  #2 

10.43 

3.96 

GaAs  #3 

10.50 

3.93 

_ Avg _ 

10.37  ±0.16 

3.98  ±  0.06 
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Conclusion 


To  conclude,  we  have  introduced  the  AATT  transverse  thermoelectrics  in  semiconductors  with 
anisotropic  electron  and  hole  dispersions,  whereby  longitudinal  electrical  currents  drive 
transverse  heat  flow,  and  deduced  the  governing  transport  equations.  The  phenomenon  is  shown 
to  be  distinctly  different  from  longitudinal  thermoelectrics  and  the  Ettingshausen  effect.  The 
optimal  transverse  figure  of  merit  ZlT  are  deduced  based  on  an  electron-hole  two-band  model. 
The  InAs/GaSb  T2SL  has  an  appropriate  AATT  band  structure  as  simulated  with  the  full-band 
envelope  function  method.  Realistic  estimates  of  the  T2SL  suggest  thermoelectric  cooling  of  8  K 
at  room  temperature,  which  is  competitive  with  recent  results  in  on-chip  cooling  by  thin-film 
longitudinal  thermoelectrics.  Estimates  of  ZlT  =  0.5  would  yield  AATT  cooling  of  56  K  at  room 
temperature,  and  thick,  exponentially  tapered  devices  could  cool  to  arbitrarily  low  temperatures 
independent  of  ZlT  value.  We  have  also  established  technological  capabilities  for  making  thin 
film  thermocouples  for  the  purpose  of  thermoelectric  calibration  of  T2SL  thin  films.  This  work 
led  to  5  published  papers,  1  preprint  (submitted  to  PRL)  and  2  patent  applications. 
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